arXiv: 1505.04446v 1 [astro-ph.HE] 17 May 2015 


Draft version May 19, 2015 

Preprint typeset using DTgX style emulateapj v. 5/2/11 


PREDICTIONS FOR REVERBERATING SPECTRAL LINE FROM A NEWLY FORMED BLACK HOLE 
ACCRETION DISK: CASE OF TIDAL DISRUPTION FLARES 

Wenda Zhang, Wenfei Yu 

Shanghai Astronomical Observatory and Key Laboratory for Research in Galaxies and Cosmology, 

Chinese Academy of Sciences, 80 Nandan Road, Shanghai, 200030, China. E-mail: wenfei@shao.ac.cn 


Vladimir Karas, Michal Dovciak 

Astronomical Institute, Academy of Sciences of the Czech Republic, 

Bocnf II 1401, CZ-14100 Prague, Czech Republic. E-mail: vladimir.karas@cuni.cz 
Draft version May 19, 2015 

ABSTRACT 

Tidal Disruption Events (TDEs) can be perfect probes of dormant SMBHs in normal galaxies. 
During the rising phase, the accretion luminosity can increase by orders of magnitude in several weeks 
and the emergent ionizing radiation illuminates the fresh accretion flow. In this paper, we simulated 
the evolution of the expected spectral line profile of iron due to such a flare by using a ray-tracing code 
with effects of general relativity (GR) taken into account. We found that the time-dependent profile 
changes significantly with black hole spin, inclination angle with respect to the black-hole equatorial 
plane, and the expansion velocity of the ionization front. At low values of spin, a “loop” feature 
appears in the line profile vs. time plot when the inclination is no less than 30° and the expansion 
velocity v exp is no less than half speed of light, due to a shadow in the emission of the truncated 
disk. In the light curve two peaks occur depending on the inclination angle. At large u exp , a shallow 
“nose” feature may develop ahead of the loop; its duration depends on the expansion velocity and the 
inclination angle. We explore the entire interval of black hole spin parameter ranging from extreme 
prograde to extreme retrograde rotation, — 1 < a < 1. In the prograde case, a low-energy tail appears 
to be more pronounced in the evolving centroid energy of the line. Our results demonstrate the 
importance to search for X-ray spectral lines in the early phase of TDE flares in order to constrain 
black hole mass and spin, as well as properties of the innermost accretion flow. 

Subject headings: accretion, accretion disks — black hole physics — galaxies: nuclei 


1. INTRODUCTION 

X-ray spectroscopy reveals that atoms and ions of iron 
can produce a strong Fe K spectral line, which is thought 
to be shaped by relativistic effects as it emerges from the 
inner regions of an accretion disk. Hard X-ray photons 
from a hot black hole corona can irradiate the accretion 
flow. Emission from an irr adiated flow then con tains the 
reflection component (see Fab ian fe Ru ss 2010, and fur¬ 
ther references cited therein for a review). The emerg¬ 
ing signal includes both the continuum flux and spectral 
lines, namely, the prominent fluorescent emission. The 
latter is produced by de-excitation of an atom with a 
vacancy in the K-shell being filled by an outer electron. 
The yield of the fluorescent line scales as Z 4 (quartic 
power of the atomic number), hence the iron fluorescent 
line is particularly strong. About 80% of the Fe K line 
photons come directly from the Fe atoms/ions without 
experiencing scattering. 

X-ray spectra have been successfully used to study 
black holes that are surrounded by an accretion disk. 
Starving black holes are more difficult to explore. Also 
when the surrounding material is inactive (e.g. a non¬ 
accreting clumpy torus) then the electromagnetic signa¬ 
tures cannot be employed to reveal the black hole. In 
such circumstances the immediate vicinity of the black 
hole horizon can be probed when a star passes too close 
to the black hole, so that it be comes destroyed in a Tidal 
Disruption Event (TDE; see iKomossa fc Greinerl 119991: 


iGezari et al.l 120091: iKocsis fe Loebl 120141 lKomossall2015L 
and further references cited therein). The resulting tran¬ 
sient accretion is delayed and it occurs much later after 
the moment of disruption of the star, as the debris mate¬ 
rial proceeds to the black hole on the viscous time-scale. 
Eventually, the gas approaches the innermost stable cir¬ 
cular orbit and at this stage the enhanced accretion gives 
us an opportunity to examine the black hole by lighting 
up its environment and producing distinct spectral sig¬ 
natures. 

While predictions for optical spectra have bee n dis¬ 
cussed in more detail (|Strubbe fc Qua taert 2009, 2011; 
iKomossal |2012[) , X-rays reach closer to the black hole 
and relativistically skewed spectral line can strongly con¬ 
strain the paramet ers, thus enab ling us to e xplore effects 
of strong gravity (iChen fc Bogdanovid l 2015 ). 

As the relativistic line is produced near the horizon, the 
line energy will be influe nced by the gravitat i onal redshif t 
and the Doppler effect (jFabian et al.l 119891: lLao a EMU). 
Hence, the observed spectral profile will be different from 
the rest-frame shape of the spectral feature, in particu¬ 
lar, the line centroid energy will be shifted (the change 
towards higher or lowe r energy is possible depending on 
parameters: lKarasll2006f h As a result, the reflected spec¬ 
tra are distorted and the observed iron Ka line profile 
becomes broad and skewed. Assuming that the flow ex¬ 
tends d own to the innermost stable circular orbit (ISCO; 
iBardeen et al.|[l972l ). the emerging spectrum can be used 
to infer the spin of the central compact object. This has 
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been done for both Gala ctic black hole X-ray binaries 
(XRBs; e.g.. iMiller et alJ 12013 ) and supermassive black 
holes (SMBHs) in Active Galactic Nuclei (AGNs; e.g., 
iRisaliti et al.ll2013l ). Besides that, the iron line profile 
can h elp us p robe the geometry of the hard X-ray emit¬ 
ter ([Fabian et alJl2QQ2h . 

The emergent spectral line depends on the ionization 
state of t h e gaseous material. As first recognized by 
I Ross fe Fabianl (j 19931 ). in the case of X-ray irradiation 
the strength of the resulting line emission is suppressed 
due to Auger destruction when the ionization parameter 
£ (defined as £ = L/nr 2 , where L is the hard X-ray lumi¬ 
nosity, n is the electron number density of the irradiated 
disk, and r is the radius) is in the range of 100-500. Then 
the equivalent width of the line increases for £ > 500 due 
to the presence of Fe XXIV and Fe XXV (the emission 
lines of iron are dominated by contributions from these 
two species). At very high ionization state (£ > 5000), 
the emission line would be suppressed again sin ce a l arge 
fract i on of the iron atoms ar e fully ionized ([Matt et al.l 
H993L [19961: Ikaras et al.ll2000l) . 

A detailed shape of the observed (relativistically 
skewed) line can in principle be used to constrain pa¬ 
rameters of the accretin g system , _in part icular, the 
black hole spin. However. [Reynolds fe Begelman ( 199 7) 
showed that there is a degeneracy in determ ining the sp in 
by time-averaged X-ray spectra. iDovciak et al.l (2 004 b) 
demonstrated that the interplay between radius of the 
inner edge (offset from ISCO) and spin can be disentan¬ 
gled, but only with spectra of very high quality. 

On the other hand, in a time-resolved process, the vari¬ 
ability of the central hard X-r ay sources would “echo ” on 
the irradiated accretion fl ow. iRevnolds et al.l ([ 19991) and 
lYoung fe Reynolds ( 2000 ) studied the resulting reverber¬ 
ation features of black hole spin, given the geometry of 
the X-ray emitting source, and they simulated the result¬ 
ing response to a S function flare in AGN. Recently, the 
promising reverberation technique has become debated 
in th e context of time lags em erging in the new X-ray 
data ([Zoghbi et al.ll2010l I2013D . 

For AGN, the lamp-post model of the irradiated in¬ 
ner accretion d isk has been adapted as a suitable work¬ 
ing hypothesis dMatt et al.lll99l( iMartocchia et al .1120001: 
iMiniutti fe Fabianl 120041 ). In this scenario the accretion 

disk is illuminated by a primary source which is located 
above the disk plane. The frequently-used assumption 
places the primary source on the common axis of the 
rotating black hole and the equatorial accretion disk. 
The source of primary photons represents a corona or 
a base of the jet. Following an onset of the illumina¬ 
tion, the primary photons travel along null geodesics un¬ 
til some of them hit the surface of the accretion disk, 
where they are reprocessed and the emerging signal (in¬ 
cluding the reflection spectral line photons) then pro¬ 
ceeds toward a distant observer. In such a geometrical 
set-up, the illuminating light front can effectively move 
at a superluminal speed across the disk surface. On 
the other hand, formation of the irradiation source in 
the lamp-post model could be enhanced by relic inac¬ 
tive (collision-less) plasma structures persisting as r em¬ 
nants of previous accretion episodes jCremaschini et al.l 
120131 : iCremaschini fe Stnchlfk 1120141: iKovar et al.ll2Q14fi T 

including the potential disruption of asteroids and small 
bodies by the supermassive black holes ([Kostic et al.l 


120091 . 120H [Zubovas et al.ll20H) . 

In the present paper, we investigate a possibility of ion¬ 
izing the freshly accreted material of debris flow, formed 
from a disrupted or partially damaged star as well as 
pre-existent gas near SMBH. This material is activated 
by X-ray illumination that would be produced near the 
black hole from a hot corona or a hot inner accretion flow 
as the accretion rate grows above the Eddington limit. 

An expanding ionization region develops and proceeds 
through the accretion flow. This is thought to be trig¬ 
gered by a TDE in a normal galaxy center, where the 
induced variation of the spectrum can be revealed more 
easily than in a typical AGN (the latter would be in¬ 
trinsically highly variable). In consequence of the chang¬ 
ing ionization the spectral line becomes modulated in 
strength and energy via increasing ionization. Indeed, 
it occurs that the medium near the central black hole is 
often in the form of a highly inhomogeneous and clumpy 
cloudlets which stay undetected unless illuminated by 
an external source. The actual flare is expected to occur 
somewhat delayed with respect to the pericenter passage 
of the affected star, as the illuminating corona builds 
gradually during the (super)-Eddington phase of fast ac¬ 
cretion following the moment of TDE. 

A non-negligible radial infall velocity of the accreted 
material is likely. The ionization front therefore moves 
radially outward at v = f exp (r) not exceeding the speed 
of light. We set t> exp — const as one of the free model 
parameters that could be constrained by the observed 
spectra. The scenario of ionization wave being induced 
by TDE and launched from the inner disk provides an 
alternative scheme with a small number of free param¬ 
eters and somewhat different properties with respect to 
changing spectral properties of the evolving line. 


2. PROSPECTS FOR IRON LINE REVERBERATION FROM 
TIDAL DISRUPTION EVENTS 


2.1. TDEs near SMBHs 

Tidal disruption events (TDEs) can be used to probe 
the cosmological SMBHs in the centers of inactive galax¬ 
ies. A TDE occurs if the stellar orbit falls within radius 
r = R p such that the ratio between stellar surface gravity 
and tidal acceleration at pericenter, 


Rl 


1/2 


< 1 , 


(1) 


where R p , R*, M* and M m are the pericenter ra¬ 
dius, stellar radius, stellar mass and the central SMBH 
mass, respectively ([Luminet & Marck 1985|; lReesl fl988: 
lEvans fc Kochanekl 119891 : iPhinnevI 1 19891 ). The exact 
value of r] for the disruption to happen depends on the in¬ 
ternal st ructure, compressibility and rotation of the star 
(lEgglctoni 1983b Simulations confirm that, after the dis¬ 
ruption, nearly half of the stellar debris should remain on 
the bound orbit and they would be eventually accreted 
by the SMBH. Based on the viscous mechanism, the mass 
fall-back rate is estimated to be 


M fb 


1M* 
3 tf b 


-5/3 


( 2 ) 


which can significantly exceed the Eddington accretio n 
rate for a period of weeks to years (iStrubbe fc Quataertl 
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l200fl IMTTT The mass of the SMBH s associated wi th 
TDE s is peaked at ~ 10 6,5 solar mases (Sto ne fc Metzgerl 
120141 ). The presence of a stellar ring or a secondary 
intermediate-mass black hole near a SMBH can enhance 
the TDE rate by increasing t he eccen t ricit y fluctua¬ 
tions of stellar trajectories (jVpkrouhlicky & Karas 1998; 
iKaras fe Subrll20l2t IKaras et al.112014]: iLi et alJl2015l ). 

In the X-ray and Optical/UV bands , ~ 20 can¬ 
didat e TDEs ha v e be e n reported ()R,enzini et al.l 
19951: iBade et all 119961: iKomossa fe Greinerf" 19991: 
Greiner et al.l l2QQ0t iDonlev et al.l 120 02: Ge zari et al l 
20061 [2008. 2009, [20H iMaksvm etTaT! I 20 I 0 I ). These 


events are characterized by thermal emission with 
temperature of ^ 10 4 -10 5 K, and the peak bolometric 
luminosity of about 10 43 -10 45 erg s _1 . For the events 
with good coverage during the decay, the luminosity 
decline was consistent with a power-law delay with 
index of —5/3. Compared with quiescent supermassive 
black holes (SMBHs), such as Sgr A* in our Galaxy 
radiating at the level of 10 35 erg s -1 , the luminosity can 
increase significantly, by several (even up to ten) orders 
of magnitude during TDE outbursts. 

In 2011, two putative TDEs were detected: Swift 
J1644+57 (e.g. , [Burrows et al.l 1201 it iLevan et~akl 120111 : 
Zanderer et al.l 1201 lh and Swift J2058.4+0516 (e.g., 
enko et aLll2012f ) . Their X-ray spectra showed a power- 
law shape, and the peak luminosity of ~ 10 48 erg s _1 far 
exceeded the Eddington limit of a SMBH of a presumed 
reasonable mass. They also exhibited radio emission with 
isotropic luminosity of ^ 10 42 erg s -1 . These are believed 
to be due to a relativistic jet viewed at a very small an¬ 
gle. The amplitude was very large in these two cases. 
For Swift J1644+57 the ROSAT 3cr upper limit was 
2.4 x 10“ 13 erg s -1 cm -2 , and the peak flux during the 


erg s 1 cm 2 (iBv 


et al.l 


erg s 

outburst was ^ 5.0 x 10 -9 
EllD), showing an amplitude of more then four orders 
of magnitude. For Swift J2058.4+0516, the ROSAT up¬ 
per limit was ^ 10“ 13 erg s -1 cm -2 while the peak flux 
during outburst was ~ 10“ 10 erg s -1 cm -2 , showing fluc¬ 
tuations with amplitude of three orders of magnitude. 

TDEs can be revealed with a better chance in sources 
where the central sup ermassive b lack ho lejs not fed by 
vigorous accretion (jSaxton fe Komossal |2012 ). In this 
context, AGN have been an unfavored type of objects 
to reveal tidal disruptions because accretion is an inher¬ 
ently variable process which can hide signatures of TDEs. 
The X-ray light curve of a typical AGN varies on mul¬ 
tiple time-scales, and so a sporadic hare by a TDE can 
hardly be discovered in the background of stochastic fluc¬ 
tuations due to the intrinsic properties of the accretion 
how. Relatively quiet low-luminosity nuclei appear to be 
more suitable for this task. 

Follow-up optical spectroscopy observations can give 
additional inf ormation on characte r istics to TDE host 
galaxies f e.g.. IKo mossa et al. 2004; S axton et al.l 120121: 
I Clausen et al . 2012; Wang et al. 20121). Here we propose 
that the evolving X-ray spectral line prohle expected in 
the very early phase of TDE accretion hares can provide 
a unique probe of the SMBHs in an otherwise inactive 
galaxy. Well after the tidal disruption happens, an ac¬ 
cretion disk forms from the circularized how of debris. 
The subsequent evolution is crucial for understandin g the 
light curves in TDE events (jShen fe Mat znerlfe 0141 ). 


The accretion of the debris material starts from low 
accretion rates and it continues ah the way to super- 
Eddington rates during the rising phase of TDEs. In 
the early rising phase, the accretion process is expected 
to proceed in the regimes of radiatively inefficient/hot 
accretion. The resulting signal mimics the one seen in 
many AGN where the lamp-post model has been applied 
to model the irradiation of the disk. A newly formed ac¬ 
cretion disk will become ionized when the source reaches 
and exceeds the Eddington luminosity. 



Figure 1 . A schematic representation of a “partially ionized ac¬ 
cretion flow”, where an ionisation wave front is launched from 
around the inner edge and proceeds outwards. The stellar disrup¬ 
tion causes enhanced accretion onto SMBH which can reach and 
exceed the super-Eddington rate. At time t after the onset of TDE 
flare, the ionised region expands to radius r — n eX p t. The inner 
disk is fully ionized and no iron line emission can emerge; only the 
outer part can produce the line emission. The fully ionized inner 
disk and partially ionized outer ring are plotted in dark and light 
gray colors, respectively. 


The idea of this investigation is based on the fact that 
an instantaneous illumination and the resulting release of 
radiation energy by the event occurring near the center 
(at the tidal radius R p of the supermassive black hole) 
shah induce an ionizing front that propagates outwards 
and changes the ionization state of the accretion disk. 
Consequently, it also influences the spectral line emissiv- 
ity in a very specific manner. The proposed approach of¬ 
fers an interesting opportunity, nevertheless, it is a chal¬ 
lenge because the iron line flux is probably weak in TDEs, 
even in underluminous nuclei with a low accretion rate 
and no prior accretion disk. 

We constructed a simplified scheme to simulate the re¬ 
sponse of the fluorescent iron line to the central TDE 
flare, where the rising edge of the TDE flare can be 
the perfect illuminating signal for reverberation mapping 
purposes. The advantage of our description is the pa¬ 
rameterization of the model by only a small number of 
free parameters. The main purpose of this work is to 
identify characteristic patterns in the line emission, not 
only limited to iron lines, that may be used to probe 
SMBH mass and spin via general relativistic effects in 
TDEs. These are thought to be rare events by their na¬ 
ture, however, the main advantage is that they can be 
used to explore SMBH in otherwise non-active galaxies 
and therefore these systems have a potential advantage 
of being more clean than AGN. On the other hand, differ¬ 
ent but related aspects of complex iron-line diagnostics 
have been dis cussed in the context o f binary black-hole 
mergers (e.g., iMcKernan et al.l 1201 3). 
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2.2. Signatures of Expanding Ionization Front 

In the following we denote the gravitational constant 
G, the mass of the supermassive black hole M, and the 
speed of light c (units of radius, time and velocity are 
GM/c 2 , GM/c 3 and c, respectively). 

Because in the rising phase of TDE accretion flares 
the continuum flux can vary in a large range, we adopt 
a corr espondingly wide range of ionization into con sider- 
ation (|Ross fe Fabian|[r993TlGoosmann et aLll2006[) . For 
a 10 6 M© SMBH, th e peak mass fail- back rate can reach 
1.5 M© yr 1 (lEvans fe Kochaneklll989[ ). correspond¬ 
ing to the luminosity of 8.5 x 10 45 erg s -1 with the radi¬ 
ation efficiency of 0.1. 

We estimate the number density of electrons using 
the radiation pr essure d ominat ed hmit of the stan¬ 
dard disk dShakura fc Sunvaevl 1119731 "): n = 1.17 x 
10 8 r 3 / 2 cm -3 . The ionization parameter £ = L/nr 2 = 
3.3 x 10 15 r -7 / 2 erg cm s -1 . Even for r = 1000, £ is still as 
large as 1.05 x 10 5 . As £ increases inward, the accretion 
flow becomes fully ionized when irradiated by the peak 
luminosity of the outburst. 

For a particular ring centered at radius r in the 
equatorial plane, the ionization parameter at time t is 
£(£) = L(t — St) /(nr 2 ), where L(t') denotes the lumi¬ 
nosity of the irradiating emission at time t'. The fluo¬ 
rescent emission of the ring is expected to vary corre¬ 
sponding to the change of the central engine, but with a 
delay St which accounts for the photon traveling from 
the center to the ring. The reverberating fluorescent 
emission emerges at t = to + St where to is the time 
of the outburst onset, and it disappears after t\ when 
£(ti) = L(ti — St)/(nr 2 ) ~ 5000. Looking at the entire 
disk, during the outburst rise, the iron line emission is 
localized to a ring composed of partially ionized medium. 
Outside the ring the photons have not arrived yet; inside 
the ring the accretion flow is fully ionized. This partially 
ionized region expands with time (Figure [l] illustrates 
this scheme). 

The expansion velocity ^ exp of the “partially ionized 
ring” depends on the height h of the illuminating source: 
f ex p = y/l + h 2 /r 2 . If the source lies near the equatorial 
plane, ?; ex p is close to speed of light (for an illuminating 
source above the equatorial plane, the expansion would 
appear even super luminal). The infall velocity cannot 
be neglected at high accretion rates when f exp would be 
equal 1 — ty, and hence less than unity. 

3. NUMERICAL APPROACH AND SIMULATIONS 
3.1. Assumptions and the Model Set-Up 

In the local frame the line-emitting region is confined 
in a ring area of the disk. In the simulation we simplified 
the region to be a ring of narrow width (about unity in 
terms of gravitational radii), and the irradiating source 
flux evolves as a step-function profile. In this scenario the 
width of the ionization ring relies on the shape of the ris¬ 
ing edge. Then we computed the emerging time-resolved 
iron line profiles, as observed by a distant observer. 

For simplicity we assumed a rotating planar disk rep¬ 
resenting a geometrically-thin Keplerian accretion disk 
formed from the TDE debris material, extending down to 
IS CO when we see the X-ray emission associated with the 
TDE. The outer boundary of the disk in our code was set 


to be 1000, which is large enough to investigate the cur¬ 
rent reverberation problem for our purposes. The radial 
power-law emissivit y index was set to 3 ( the Newtonian 
value at large radii; I Vaughan et all l2004j h We assumed 
the iron line emission to be locally isotropic. Detailed 
analysis of the angular distribution of t he line emission 
(jSvoboda et al.l 120091: iGarcia et al.ll2014l ) showed that it 
is likely to exhibit a weakly limb-brightening dependence, 
with maximum ratio of only ^2 compared to an isotropic 
assumption. This is less extreme then I oc ln(l + l//i), 
which is another commonly assumed law for t he angula r 
dependence of an X-ray illuminated slab (lHaardt Ill993h . 
Hence, our assumption of isotropic emission would not 
change the result very much. We assumed the line to be 
intrinsically narrow and its profile to be a ^-function in 
energy centered at 6.4 keV, i.e., the laboratory energy 
of Fe I Ka. The results for spectral lines at different 
energies can be inferred from our simulation. 

With all the assumptions above, the local emissivity is 

const * r _3 ,n ex p *t-l <r < n exp *t 

\0 ) 

0, otherwise. 

This is obviously a rather simplistic parameterisation of 
a much more complicated situation, nonetheless, we find 
it useful for the sake of demonstrating the role of GR ef¬ 
fects in the context of wave perturbation that is launched 
from the inner rim of the disk and then starts travelling 
outwards. 

3.2. Ray tracing of the emission from an expanding 
ionization region near a black hole 

The basic properties of black hole X-ray spectra from 
the innermost regions of an accretion disc are w ell-known 
(jFabian et al.l l2QQ0L iRevnolds fe Nowakll2QQ3h . Gravita¬ 
tional effects act on the spectral features from black-hole 
accretion disks and coronae by smearing their spectral 
features and moving them across energy bins. In this 
way gravity exerts the influe nce_o n the_ultimate_form of 
the observed spectrum ([C unnin gham fe Bardeenl 119731 : 
iCunninghaml 119751: [K arasll200if b The reprocessed radi¬ 
ation reaches the observer from different regions of the 
system. Furthermore, as strong-gravity plays a crucial 
role, photons emitted at the same moment may even fol¬ 
low multiple separate paths, joining each other at the 
observer at different moments. Individual rays experi¬ 
ence diverse time lags for purely geometrical reasons, for 
relativistic time-dilation, and by repressing within the 
medium. 

To simulate the evolution of the observed spectral line, 
we have employe d the general- r elativis t ic ray- tracing 
method (KYcode; iDovciak et al.l l2004ai l2QQ4bl ). This 
code uses pre-calculated information about the form of 
light rays (i.e., the geometrical shape and the timing re¬ 
lations along null geodesics, which describe the photon 
propagation in the approximation of geometrical optics). 
A fine grid is used with multiple scales that allow us 
to reproduce the effects of varying energy shifts, grav¬ 
itational lensing near caustics, frame-dragging near the 
ergosphere, and time delays in the gravitational field of a 
rotating (Kerr) black hole (dimensionless spin parameter 
-1 < a < 1). 

First, the source of the emerging signal is represented 
by an expanding ionization wave launched in the equato- 
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Figure 2. Contour lines of redshift factor g{r , <j>) = F? 0 bs/-^loc in the equatorial plane of a rotating black hole for several values of its spin 
(rows) and of observer’s inclination (columns). In this plot the black hole rotates in both clockwise (a > 0) and counter-clockwise (a < 0) 
sense, as indicated (a = 0 is the static case with zero angular momentum). The radiating atoms follow Keplerian circular motion around 
the black hole, while a distant observer has the line-of-sight directed downward. In each panel the black circle in the middle indicates the 
location of ISCO. The value of g is color coded, as indicated by the color bar to the right of the figure. 


rial plane, as introduced in the previous section and fur¬ 
ther detailed below in Sec. HJ The resulting spectra and 
the corresponding light curves are produced in the form 
predicted for a distant observer, i.e., as they are expected 
to be recorded at the detector plane (defined by the view- 
angle inclination with respect to the black hole rotation 
axis) at radial infinity. The code uses the relevant func¬ 
tions stored in FITS files, so that the subsequent com¬ 
putations of the spectral changes are performed in very 
efficient manner. This allows us to reconstruct the im¬ 
pact of energy shifts (both the Doppler red/blue shifts 
due to positive/negative component of motion along the 
line of sight, as well as the overall gravitational red shift 
near the horizon), light focusing, and the photon travel 
time from the point of emission towards an observer far 
from the source. 

We have set-up a grid of the interpolation mesh 
with respect to relevant variables: black-hole spin, the 
inclination angle, and the expansion velocity of the 
front. For the spin we select exemplary values of a = 
—0.998, —0.5, 0, 0.25, 0.5, 0.75, and 0.998. Both prograde 
(co-rotating, a > 0) and retrograde (counter-rotating, 
a < 0) cases are investigated in this work. This covers 
the entire range of spin values that allow for the black 
hole horizon, while assuming that the inner rim is linked 
with the ISCO radius (the latter recedes in the retro¬ 
grade case, and so the relativistic effects become less im¬ 
portant) 0 For the inclination angle we adopt the grid 

1 In principle, the values of \a\ > 1 should not be excluded from 


values in the range from i = 0 (corresponding to the ro¬ 
tation axis) up to i = 85° (almost edge-on view), with 
the resolution of 5°. For the expansion velocity, we set 
the reference value at v exp = 1 (corresponding to the 
case of an illuminating source lying in the disk plane and 
the ionization taking effect immediately), and we show 
r?exp = 0.5 and 0.75. These allow us to take the effects of 
accreted debris radial inflow into account. For reference 
and comparison we also calculated the case of v exp = 1.2 
and 2.0 which can represent the super luminal expansion 
for an illuminating source located above the disc plane, 
such as the case of lam p-post scenario in Kerr metric 
(|Martocchia et ahll2000f ). 

Let us note that the constant value of the expansion 
velocity ^ exp = const is introduced to illustrate the ex¬ 
pected dependencies and it serves as a simplified param- 
eterisation. In a realistic description the ionisation front 
will proceed in the interaction with the inflowing ma¬ 
terial, which exhibits a radially dependent infall veloc¬ 
ity. Also, the relativistic effects will modify the simplistic 
scheme. Nevertheless, even the lamp-post scenario finds 
its phenomenological substantiation as it can mimic the 
role of the inner accretion disk self-irradiation, where the 
light-bending causes an additional energy input imping¬ 
ing onto the disk surface predominantly over the location 
on the disk axis at distance of roughly the light circle 

the analysis a priori , however, we do not consider this possibility 
because then a naked singularity develops and this would bring 
us beyond the standard scenario of black-hole accretion (work in 
progress). 
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Figure 3. Distribution of the relative time delay of photons originating from the black hole equatorial plane. The set-up of the plot is the 
same as in Fig. [2] and the values of the time delay are encoded by the color bar (values are given in units of GM/c 3 ). Notice the prevailing 
influence of strong gravity near the horizon, where the light travel time quickly grows above the corresponding Newtonian (flat space-time) 
value. This effect distorts the perfectly straight (horizontal) of the Newtonian time-delay contours which in the relativistic case encircle 
the horizon. 


(3 GM/c 2 for a = 0). For this reason we think it is useful 
to explore the entire range of v exp from small values to 
greater than unity. 

4. RESULTS 

We show how the different values of black hole spin, 
inclination angle and expansion velocity affect the evolu¬ 
tion of the relativistic spectral line emission as predicted 
in this model. 

To better understand where the features in the time- 
resolved line profile originate, we first illustrate the dis¬ 
tribution of redshift and time delays for photons emitted 
in the equatorial disk in Kerr space-time for different val¬ 
ues of black-hole spin a and observer’s inclination angle 
i. These functions then determine the time-dependent 
profiles of the iron line. 

Figure [2] shows the distribution of the redshift on the 
equatorial disk plane, for different values of a and i. The 
redshift g(r, </>) = E 0 \y S /E\ 0C is affected by both gravita¬ 
tional field of the central black hole and also the Kep- 
lerian circular motion of the disk. The redshift function 
reaches minimum in the innermost region of the accretion 
flow, which is closest to the black hole (see the bottom 
panels of Fig. [2j where a = 0.998 and the ISCO lies on 
r = 1.23). The redshift factor can be as small as ^ 0.2 
in the inner region around ISCO. Photons emitted from 
this area contribute to the extended red wing of the broad 
iron line in X-ray spectra of AGNs and galactic XRBs. 
The Doppler effect is more profound for large inclination 
angle edge-on systems, as the component of the velocity 


along the line-of-sight increases with the inclination. 

Fig. [3] shows the effect of purely geometrical (general 
relativistic) relative time delay of photons arriving from 
the equatorial plane. In the case of prograde rapidly 
rotating black hole (a = 0.998, the bottom panels), the 
spectral-line photons experience large time-delays in the 
inner disk. For the high inclination systems (i = 85°; 
right panels), the photons from behind the black hole are 
delayed when they pass the black hole. In the extreme 
case of both large inclination and rapidly rotating black 
hole (the right-bottom panel), the frame-dragging effect 
is also seen. 

In Figures 0H6] we show the time-resolved line pro¬ 
file, both for Schwarzschild (a = 0) and extreme Kerr 
black holes (a = 0.998,-0.998). We set zero time as 
mean time, defined by to = f tf(t)dt/ f f(t)dt where 
fit) denotes the photon counts recorded at the detector 
plane 0 The corresponding light curves are shown in fig¬ 
ure [7J where Schwarzschild and Kerr cases are plotted 
in dash-dotted and solid lines, respectively. Let us note 
that here we show the (computed) light curves of the 
background-subtracted model; hence, the plot exhibits 
the expected variation of that part of the observed sig¬ 
nal restricted to the energy of the iron line emission. 

To characterize the changing profile we employ the line 

2 To compare with observations, it would be natural to define 
time zero to be the moment when the onset of the TDE flare reaches 
the observer. However in our assumption the illuminating source 
lies in the center of the disk plane which is inside the event horizon, 
so it is actually not feasible to define zero time in this way. 










































Reverberation light curves and spectral lines from a black hole accretion disk 


7 


> 

0 


cn 

0 

C 

LU 


. vP i .. 

. ?*?, .:.. 

.. p.r. . .3 . 

.,,^P.~ . mssrn 




.1 



i 

^1 

: 

: : 

....... : 


w 

; : 

/%. ; 

; ’ ./V- ; 

^ . 


-60 -40 -20 0 20 -40 -30 -20 -10 0 10 -40 -30 -20 -10 0 10 -40 -30 -20 -10 0 10 20 


Time (GM/c 3 ) 


Figure 4. Evolution of the observed (computed) iron line photon counts (background subtracted) for Schwarzschild black hole (a = 0), 
for different inclination angles and expansion velocities. The photon count values are normalized to the maximum value for all panels and 
encoded by gray-scale (in arbitrary units). The corresponding inclination angle and expansion velocity are shown at left and top axes, 
respectively. 


centroid energy, 


E, 


f E f(E) dE 
f f(E) dE 


(4) 


where E is the photon energy, and f(E) is the photon 
number density flux in the line, both of which are mea¬ 
sured with respect to the observer’s frame. While the 
details of the evolving spectrum are telling about GR 
effects from the theoretical point of view, the centroid 
energy provides a robust characteristic that can be com¬ 
pared with observation. 

For the Schwarzschild black hole case at i = 5° (i.e. 
viewed almost along the symmetry axis), the iron line 
has a narrow profile at any time and the line energy in¬ 
creases with time asymptotically towards the rest-frame 
value. As the inclination angle increases from 5° to 85° 
(close to an edge-on view), the observed line profile at 
any time becomes broadened due to an increasing veloc¬ 
ity component of the disk material along the line of sight. 
Once f ex p >0.5 and i > 30°, a loop is seen in the line 
profile vs. time plot (see Figs. 0H6]). The time span of the 
loop increases with the inclination angle. At ?; ex p > 1.0, 


two peaks can be seen in the line light curve, as shown in 
Fig. [71 and the first and the second peaks correspond to 
the left and the right edges of the loop, respectively. 

We computed the time span of the two peaks and we 
found that their time difference stays roughly constant 
with varying expansion velocity. For example, when i = 
60°, the time separation is 13.8, 13.6, 13.4 for f exp = 1.0, 
1.2 and 2.0, respectively. On the other hand, the time 
separation and intensity ratio of the two peaks depend 
on the inclination. The mean time separation = 

1.0, 1.2 and 2.0) is 7.1, 13.6 and 17.3 for i = 30°, 60°, 
and 85°. The intensity ratio of the first to the second 
peak decreases with the inclination angle in most cases. 
During the first peak, a dip occurs in the centroid-energy 
evolution, and between the two peaks the centroid energy 
exhibits a bump (see Figure [8)). 

The time delay between the central TDE X-ray irra¬ 
diating flare and the induced iron line is denoted as At 
and contributed by two factors: the first one is the ex¬ 
pansion time, £ exp ~ r/v exp , and the second one is the 
travel time needed from the disk to the observer, tdelay 
(Fig. [3]). At large radius where the spacetime is flat to a 


























Wenda Zhang et al. 


> 

0 


cn 

0 

C 

LU 



Time (GM/c 3 ) 


Figure 5. As in the previous figure, but for a fast-rotating black hole (spin a = 0.998). Prograde sense of rotation is assumed, i.e., the 
accreted material co-rotates with the black hole. 


good precision, tdeiay — — rsin</> + const, where </> is the 
azimuthal angle on the disk plane. As t; exp is of the order 
of unity, the two are comparable, and At is dominated 
by the longer one. Hence in the case of t; ex p > 1, At is 
dominated by tdeiay- 

At large inclination angles, the photons originating 
from the disk in front of the black hole reach the ob¬ 
server first, then come those photons from the disk near 
ISCO, and the last are those photons from the disk be¬ 
hind the black hole. The loop feature in the line profile 
vs. time plot is seen due to the existence of a hole (due to 
our assumption that the disk would be truncated inside 
the ISCO). The left edge of the loop corresponds to the 
innermost isochronous disk plane unaffected by ISCO in 
front of the black hole, and the right edge of the loop cor¬ 
responds to the innermost isochronous plane unaffected 
by ISCO behind the black hole. Roughly, the time span 
of the loop equals to l ~ 2risco sin i (where nsco is the 
radius of the ISCO), which is consistent with the time 
separation between the two observed line intensity peaks 
in the light curve. 

The relative intensity of the second peak is enhanced 
for large inclination angles as it is formed mainly by the 
photons emerging from the disk behind the black hole 


(that get more focused at high inclinations). In the case 
of low expansion velocity, At is dominated by £ exp , and 
the observed image of the emission region expands out¬ 
wards from the ISCO. In the cases of high expansion 
velocity combined with high inclination angle, there are 
“noses” occurring ahead of the loop. In these cases, the 
observed image of the emission region first expands out¬ 
side the ISCO and only later on reaches the region below 
the ISCO. This expansion creates the nose that is fol¬ 
lowed by the loop and two intensity peaks are visible in 
the light curve at the time when the image reaches and 
subsequently leaves the ISCO. 

The duration of the nose increases with inclination an¬ 
gle, and with the expansion velocity as well. The noses 
are contributed by photons emerging from the disk in 
front of the black hole. Higher inclination and larger ex¬ 
pansion velocity cause that the first image of the emission 
region comes from further away from the centre which 
prolongs the nose in time-energy graphs. The centroid 
energy of the noses is close to the rest frame energy of 
6.4 keV (this is also reflected in Figs. 0H6|). 

Comparing with the Schwarzschild case, in extreme 
Kerr black hole cases, the line energy can extend to both 
lower and higher energies, due to stronger gravitational 
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Figure 6. As in the previous figure, but for a counter-rotating black hole (spin a = —0.998). Unlike the previous figure, retrograde 
rotation means that the accreted material has its angular momentum opposite to the black hole. This leads to receding ISCO radius, and 
therefore decreasing the span of energy that the event covers in the plot above in comparison with the case of prograde rotation. 


and Doppler energy shifts as the inner edge shrinks with 
increasing spin. The loop feature is no longer seen, while 
the nose stays unaffected. 

At v exp >0.5 speed of light, a low-energy tail appears 
in the profile (see Fig. [5j). A similar signature ha s been 
noticed for Kerr black holes (lYonng fc Reynolds 2000 ). 
As the inner edge is very close to the black hole, the pho¬ 
tons from this region of the disk suffer from large time 
delays (see bottom panels of Fig. [3]) and strong gravita¬ 
tional redshift (see bottom panels of Fig. [2]), and for this 
region At is dominated by tdelay at v ex p > 0.5. These 
two contributing effects produce the delayed low-energy 
tail. As a result, the two intensity peaks do not develop 
in the light curve, while a dip remains to be seen in the 
centroid energy evolution. 

5. DISCUSSION AND CONCLUSIONS 

We simulated the time-resolved iron-line profiles from 
TDE flares, as they are expected to be seen in a distant 
observer’s frame. Because the maximum mass fall-back 
rate can be supper-Eddington in TDEs, a large range 
of luminosity is expected in these events. Hence, the 
energetic photons from the central region can bring the 


accretion how to high ionization state. The rising edge 
of TDE hares can produce and modulate the fluorescent 
iron line. The line emission would disappear in the ac¬ 
cretion how where the disk material is highly ionized (in 
terms of £ > 5000). 

We explored a simple parameterization according to 
which the velocity of a moving ionization front t> exp de¬ 
fines a different kind of variable compared to the more 
standard height h of the lamp-post scenario. Other pa¬ 
rameters of the model are the black hole spin a (Kerr 
metric is assumed) and the line-of-sight view angle i 
of the system. Motivation of the model set-up arises 
in the context of spectral-line reverberation mapping of 
the spectral line from the debris disk illuminated after 
a tidal-disruption event. We suggest that this scheme 
can be useful to study the properties of tentative TDEs, 
where the emerging spectral line is modified by an out¬ 
ward propagating circle which grows and changes the line 
emissivity. The accretion process is expected to be highly 
non-stationary in such a situation. 

We conclude that notable signatures can be expected: 

(i)—For i > 30°, f exp A) 1? and small or negative (ret¬ 
rograde) spin, a loop feature appears in iron line profile 
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Figure 7. The observed count rates in the fluorescent iron line (as computed by our ray-tracing method, background subtracted), for 
different inclination angles i and expansion velocity u eX p, as indicated at each corresponding row and column, respectively. The dash-dotted, 
solid, and dotted lines correspond to the extreme retrograde Kerr (a = —0.998), non-rotating Schwarzschild (a = 0), and extreme prograde 
Kerr (a = 0.998) black holes, respectively. Each lightcurve is normalized by its maximum value. 


vs. time plot and it manifests itself in the line flux show¬ 
ing two peaks in the corresponding light curve. The time 
separation and intensity ratio of the two peaks mainly 
depend on the inclination angle. For large values of pos¬ 
itive (prograde) spin parameter, the radiation emitted 
close above the horizon destroys the loop (see also the 
point (iii) below). 

(ii) — For large expansion velocities of the ionization 
circle and large inclination angles, a nose appears ahead 
of the loop in the iron line profile vs. time plot. The 
duration of the nose is determined by both the expansion 
velocity and the inclination angle. 

(iii) — For v exp > 0.5 and a 1 (extreme prograde 
rotation), a more prominent low-energy tail is seen. Oc¬ 
currence of this feature can serve as a probe of the fast- 
rotating Kerr black hole. On the other hand, in the non¬ 
rotating and retrograde cases the low-energy tail is sup¬ 
pressed and it happens earlier in time. 

Let us note that the retrograde sense of rotation is 
often dismissed in accretion disk modelling and spectra 
fitting because it appears to be less likely compared to 
the prograde case. It can be argued that the long-term 
evolution the accretion flow tends to align the angular 
momentum vector of the flow with that of the black 
h ole. However, t he si tuation with TDEs is quite different 
(|Sochora et al.lfJoilD . An infalling star approaches the 


central SMBH from the surrounding nuclear star-cluster, 
which can be considered as spherical and isotropic in 
zero approximation, and so counter-rotating motion is 
equally possible as co-rotating with the black hole. Fur¬ 
thermore, counter-rotating stellar rings and disks have 
been reported in some galaxies, so we included this pos¬ 
sibility in our considerations. 

Once a TDE-induced relativistic line is detected and 
its parameters reliably determined in the future, these 
signatures in time-resolved signal can be used to con¬ 
strain the black hole spin, inclination, and the expansion 
velocity of the ionization front. 
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Figure 8. The expected time variations of iron line centroid energy E cen for different inclination angles i (rows) and expansion velocities 
r exp (columns). The line types as in the previous figure: the dash-dotted line for a = —0.998), the solid line for a = 0, and the dotted 
line for a = 0.998 cases, respectively. The corresponding inclination angles i and expansion velocities i; eX p are also given. The rest-frame 
energy 6.4 keV of Fe I Ka spectral line is indicated by the horizontal dashed lines for the reference. 


titled “Spectral and Timing Properties of Cosmic Black 
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